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Abstract 

We reconsider thermal production of gravitinos in the early universe, adding to 
previously considered 2^2 gauge scatterings: a) production via 1 —> 2 decays, 
allowed by thermal masses: this is the main new effect; b) the effect of the top 
Yukawa coupling; c) a proper treatment of the reheating process. Our final result 
behaves physically (larger couplings give a larger rate) and is twice larger than 
previous results, implying e.g. a twice stronger constraint on the reheating temper- 
ature. Accessory results about (supersymmetric) theories at finite temperature and 
gravitino couplings might have some interest. 
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1 Introduction 



We compute the abundance of gravitinos thermally produced in the early universe at temper- 
ature T. In the usual scenario where sparticles around the weak scale keep it naturally small, 
this process implies an important constraint on the maximal reheating temperature, possibly 
saturated if such gravitinos are all observed Dark Matter (DM). If instead sparticles exist much 
above the weak scale, gravitino production is one of their very few experimental implications 
that survive. 

The gravitino production thermal rate was previously computed in [1, 2] at leading order 
in the gauge couplings (and g2,gY in [3]; we will add effects due the top Yukawa coupling, 
which also has a sizeable value). This roughly amounts to compute 2 — >■ 2 scatterings (like 
gluon + gluon — > gluon — > gluino + gravitino), with thermal effects ignored everywhere expect 
in the propagator of the virtual intermediate gluon: a massless gluon exchanged in the t-channel 
gives an infinite cross-section because it mediates a long-range Coulomb-like force; the resulting 
logarithmic divergence is cut off by the thermal mass of the gluon, m ~ gT, leaving a In T/m. 
The explicit expression for the number of scatterings per space-time volume, at leading order 
in the dominant QCD gauge coupling, was found to be [2, 3]^ 

where Mpi = 2.4 10^^ GeV is the reduced Planck mass, M3 is gluino mass and 7/13/2 is the 
gravitino mass. This production rate unphysically decreases for (73 > 0.7 becoming negative for 
g3>1.2. Fig. 1 shows that the physical value, gs ~ 0.85 at T ~ lO^'^GeV, lies in the region 
where the leading-order rate function /{g^) (dashed line) is unreliable. Fig. 1 also illustrates 
our final result (to be precisely described in section 4.2): / will be replaced by the continuous 
lines, which agree with the leading order result at g ^ m/T <^ 1 and differ at (7 ~ 1. 

Let us now explain why the leading-order approximation in (1.1) starts to be inadequate 
already at (7 ~ 0.7. In thermal field theory higher order corrections are usually suppressed by 
g/n: somewhat worse than the usual expansion coefficient (g/Tf)^ at T = 0, but still typically 
good enough at g ~ 0.7. Naive power counting fails (without signaling a breaking of the 
perturbative expansion) when some new phenomenon only starts entering at higher orders, and 
this is what happens in the case of gravitino production: a new simpler process gives corrections 
of relative order {gT^Y- The gravitino couples to two particles with different thermal masses: 
gluon/gluino, and quark/squark. Since thermal masses grow like T, this gives rise to a new 
process with a rate growing like T^: gravitino production via decays, such as gluon gluino 
-|- gravitino, whose rate can be crudely estimated as 

mT^T m'^T^ g^ 

Tdecay ~ ~ ~ (1-2) 

Indeed 7decay is of course proportional to the decay rate at rest F ~ m^/vrMpj; which is slowed 
down by the Lorentz dilatation m/T factor; the takes care of dimensions, and less vr are 
present at the denominator because a 1 ^ 2 decay involves less particles than a 2 — > 2 scatter- 
ing. So, despite being higher order in g, the decay rate can be enhanced by a phase space factor 



^Since we will adopt a different technique, we cannot resolve the minor disagreement between the results 
of [2] and [3]. Notice also that, for later convenience, in eq. (1.1) we explicitly show the power tt^ (following 
from the phase space for scattering processes, and dictated by naive dimensional analysis), which is explicitly 
present in [1] and partially hidden in numerical coefhcients in [2]. 
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Thermal vector mass, m/T 

Figure 1: Functions fs, /2 and fi that, as summarized in section 4-2, describe our result for 
the gravitino production rate from SU(3)c (upper continuous curve, in red), 811(2)^ (middle 
continuous curve, in blue), U(l)y (lower continuous curve, in green) gauge interactions. The 
arrows indicate the MSSM values of the thermal mass atT 10^ GeV. The lower dashed curve 
shows the result from [2], which agrees with our result in the limit of small gauge coupling, and 
behaves unphysically for relevant 0{1) values of the MSSM gauge couplings. 



TT^. Subsequent higher order corrections should be suppressed by the usual g/n factors. Our 
goal is including such enhanced higher order terms, and this finite-temperature computation is 
practically feasible because a decay is a simple enough process. 

So far we explained the physical picture in a simple intuitive way. A more precise technical 
language is necessary to present how we will proceed. To get the gravitino production rate 
we actually compute the imaginary part of the gravitino propagator in the thermal plasma. 
Thermal effects distort the dispersion relations E{k) of gluons, gluinos, quarks, squarks by i) 
adding a thermal mass E"^ = k"^ + m?{k) to the modes already existing at zero temperature; ii) 
by introducing new collective excitations (gluons with longitudinal polarization, gluinos with 
'wrong' helicity, ...) with their own dispersion relation; iii) beyond the two poles mentioned 
above, the spectral densities of particles in a thermal plasma also develop a 'continuum' contri- 
bution, that can be thought of as a parton-like distribution, with a continuum range of masses. 
Physically it arises because particles can exchange energy with the plasma. 

In previous works [2, 3] the gluon thermal mass was taken into account to regulate infra- 
red divergences encountered in scattering rates, and the contribution of the gluon 'continuum' 
was computed using a standard technique introduced in axion computations [4], that allows to 
extract the rate at leading order in g. This was achieved by introducing an arbitrary splitting 
scale k that obeys gT <^ k <^T. 

We will not use this technique: because its validity is doubtful for g^, ~ 0.85, and because we 
actually want to include the enhanced higher order terms, taking into account that a gravitino 
(unlike an axion) couples to two particles with different thermal masses. We will instead 
compute the decay diagram (D in fig. 2) using resummed finite-temperature propagators for 
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Figure 2: Some Feynman diagrams that contribute to the imaginary part of the gravitino prop- 
agator. Thick lines denote resummed thermal propagators for the gluon g and gluino A. We 
do not plot diagrams involving quarks q and squarks q, but they are of course included in our 
computation. 




Figure 3: Two-loop Feynman diagrams that appear in the expansion of diagram D, that resums 
all higher loop diagrams with iterated one-loop corrections to gluon and gluino propagators. 



gluons, gluinos, quarks, squarks. The perturbative expansion of this diagram D contains the 
two- loop diagrams in fig. 3: their imaginary parts correspond to well-defined combinations of 
scattering processes, as dictated by cutting rules. This fixes how scatterings must be subtracted 
in order to avoid overcountings of effects already described by thermal masses via diagram D. 
In section 2 we compute the subtracted scattering rates, in section 4 we compute the gravitino 
production rate via 'decay', and in section 5 we add the rate due to the top quark Yukawa 
coupling. 

In section 6 we sum these effects and compute the gravitino abundance writing a set of 
Boltzmann equations that describe the reheating process, previously approximated assuming a 
maximal temperature equal to the reheating temperature Trh- Our results are summarized in 
the conclusions, section 7. 

In the passing we address some issues related to finite-temperature and to supersymmetry. 
In section 3 we list explicit values for thermal masses for all particles and sparticles, noticing 
that they obey some supersymmetric relation. Appendix A gives a (non uselessly) fully precise 
summary of gravitino interactions, and in appendices B, C we collect full expressions for the 
thermal corrections to vector and fermion propagators. 
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Figure 4: Feynman diagrams that contribute to gg A\l/ scatterings. 




Figure 5: Feynman diagrams that contribute to gX g"^ scatterings. 



2 Subtracted scattering rate 

Gravitinos with momentum P = {E,p) are produced via their couphng ^ ^S'^ /2Mp\, where 
S'^ is the supercurrent of the visible sector of a supersymmetric theory, here assumed to be 
the MSSM. The visible sector is thermalized, while the gravitino is not, since its coupling to 
the MSSM plasma is weak. According to the general formalism of thermal field theory [5], 
the production rate of such a weakly interacting fermion is related to the imaginary part of its 
propagator 11 as 

'y~,--^l''PMEnmu = j,Pu<(n dP^^. (2.1) 

Here 11^ is the non time-ordered gravitino propagator summed over its polarizations i.e. traced 
with the gravitino polarization tensor 11^,^ (appendix A gives explicit expressions): 

n<(P) = i]^Tr [U,^{P){S-'iPm-P))T] (2.2) 

where (■ ■ ■ )t denotes thermal average. We employ 11^ because it gives slightly cleaner formulae 
than Imll. Eq. (2.1) is valid at leading order in the gravitino coupling Mp/, and to all orders 
in the MSSM couplings, gY.2,3 and A^. Extracting predictions from (2.1) is limited only by our 
ability to evaluate (2.2). 

Thermal field theory cutting rules allow to see that, at leading order in the MSSM couplings, 
eq. (2.1) is equivalent to summing rates for the various tree-level processes that lead to gravitino 
production. At tree level this formalism is more cumbersome than a direct computation of 
production rates. However, in this paper we want to take into account finite temperature 
corrections to the MSSM particle propagators arising at one loop level: eq. (2.1) becomes 
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more convenient because it cleanly dictates how one must resolve ambiguities encountered in 
scattering computations that arise because Lorentz invariance is broken by the thermal plasma. 

Fig. 2 shows some main Feynman diagrams that contribute to Imll. What we actually 
compute in this paper is the first one-loop diagram 'D', using the Feynman gauge resummed 
finite-temperature propagators for the gluon and gluino in the loop. Therefore it describes a 
sum of an infinite number of multi-loop diagrams: the lowest-order ones are shown in fig. 3. 
Resummation is needed because thermal effect drastically change the gluon and gluino propa- 
gators, in particular opening a phase space for decays, such as g ^ X"^ and/or A — > g"^. Clearly 
diagram D contains this decay process. However, by cutting fig. 3, one sees that diagram D 
also describes some of the 2—^2 scattering processes computed in previous analyses [2]. There- 
fore, before starting the computation, we clarify this issue showing how the total gravitino 
production rate is obtained. 

The total scattering rate is the sum of various 2 — 2 processes A, B, C,. . . , listed in table 1. 
Each process is the modulus squared of the sum of a few amplitudes, corresponding to the 
single Feynman diagrams, that we label as s, t, u, x: 

Tscattering 

This notation indicates that often 4 diagrams contribute to a given process: 3 diagrams are 
generated by s, t and w-channel exchange of some particle among two vertices {g^ and 1/Mpi), 
and a fourth diagram arises from a quartic supergravity vertex with coupling g/M-p\. Fig. 4 
and 5 show concrete examples of the 4 diagrams that contribute to gg \&A and to g\ — >■ '^g 
scatterings, respectively. The latter rate is logarithmically infra-red (IR) divergent, because 
diagram Bj is mediated by t-channel gluon exchange, that describes a Coulomb-like scattering. 

The main result can be obtained by careful visual inspection of cutting rules: diagram D 
describes the sum \As\'^ + |Ap + \Au\^ + + \Bt\'^ + + ■ ■ ■ of the modulus squared 
of all 2 — > 2 diagrams that contain the gauge coupling g^. 2^2 scattering rates generated 
by supergravity quartic vertices are instead described by diagrams like S3 in fig. 2. Some cuts 
of the two loop diagrams like Si and S2 describe the interference terms among the various 
Feynman diagrams. (Notice that the imaginary part of a single two-loop diagram can describe 
contributions to different scattering processes). Other cuts of these diagrams describe one loop 
corrections to the gravitino vertices, that do not give any leading order contribution if thermal 
masses are neglected. Thermal masses open a phase space for 1^2 processes (we can neglect 
the decay rate generated by zero-temperature masses m, since we are interested in T ^ m), 
and we will later argue that we can still neglect thermal corrections to the gravitino vertex. 

In conclusion, the total gravitino production rate due to gauge couplings will be computed 

as 

1 = 1D + lT\ (2.3) 

the sum of diagram D (that describes decay plus modulus squared of many single 2^2 
diagrams) plus the set of remaining 2 — > 2 rates, obtained by subtracting from the total 
scattering rate 75 the effects already included in ■ju. Explicit results for 7^) and for 7!^*^ will be 
given in eq. (4.6) and eq. (2.6) respectively, and the conclusions will describe how to use them. 

Before proceeding to actual computations, we have to clarify the issues of gravitino coupling 
and gravitino gauge invariance. We are interested in T ^ m, where m denotes sparticle or 
gravitino masses: gravitino Goldstino equivalence (appendix A) means that at leading order 
in m/T the massive gravitino field can be replaced with two massless field: a massless 
gravitino coupled to the supercurrent (given in eq. (A. 28), it can be evaluated in the 
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Table 1: Squared matrix elements for gravitino production in units of gj^/Mpi{l + 
M^/3my2) ! summed over all polarizations and gauge indices. The result hold for all three 
factors of the SM gauge group with N = {1,2,3}, although the notations are appropriate for 
the SU(3) case: g,X,q,q denote gluons, gluinos, quarks, squarks. The gauge factors Cn and 
C'^ are defined in the text. Rates A and B are the sum of vector and ghost contributions. 

supersymmetric limit ignoring soft terms) plus a massless Goldstino x, coupled to the divergence 
of the supercurrent (given in eq. (A. 29), only the soft terms factored out are relevant): 

^int = + X {d,S^) _ ^2.4) 

2Mpi ■\/6^pi^3/2 

The gravitino production rate is given by ^) — 7(^/i) + lix)- While the total rate is 
gauge independent (vectors have SU(3)c ® S\]{2)l ® U(l)y gauge invariance; the computation 
of jiipfj.) also involves gravitino gauge invariance), its splitting in resummed and not-resummed 
contributions is not. We are resumming a well-defined class of effects, but we cannot system- 
atically include all the effects up to a given order in g: therefore our result has a residual 
gauge-dependence, of relative order g'^/ir'^, due to partial inclusion of higher-order terms. To 
make the computation feasible, we choose for vectors the Feynman gauge, and for the massless 
gravitino ^/^^ the gauge where its propagator and polarization tensor does not involve terms 
containing or P^, eq. (A. 23): 

nli' = -^^,riu - rv,u. (2.5) 

One first motivation for this choice is that, in the supersymmetric limit, the full supercurrent 
satisfies P^S^^ = 0, while sub-sets of 5"^ are not separately conserved: with choice (2.5) we 
never have to deal with such terms. Of course, the same gauge is used for computing both the 
resummed diagram D and the subtracted scattering rates. 

Table 1 gives explicit values for the subtracted massless gravitino and Goldstino scattering 
rates due to gauge interactions. It is important to notice that, unlike the total rate, the 
subtracted rates are infra-red convergent: no 1/t factors appear because all divergent Coloumb- 
like scatterings, like |-Btp, are included in diagram D, that we compute using thermal masses 
that provide the physical cut-off. Unlike in the conventional technique [4] employed in [2, 3], our 
technique does not need to introduce an arbitrary splitting scale that satisfies the problematic 
conditions gT <^ fc* <C T in order to control infra-red divergences. Some contributions to 
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subtracted scattering rates turn out to be negative, but the total rate will be positive and 
dominated by diagram D. In Feynman gauge, rates for the processes A and B (the ones that 
involve two vectors) actually are the sum of scatterings involving two vectors (four diagrams, 
computed with the Feynman polarization tensor ^ e^e* = —rj^y) plus scatterings containing 
two ghosts (one diagram, negative 

A curious fact happens. Despite the fact that the massless gravitino ip and the Goldstino 
X have different couplings (in particular the Goldstino has no coupling to quark/squark, and 
consequently a reduced set of Feynman diagrams), the differential production cross sections 
for these two particles are the same, process by process, up to the universal factor Mj^ /Sniy^y 
where Mi 2,3 are the gaugino masses. We don't know if there is a simple generic reason behind 
this equality. The second reason for choosing the gravitino projector of eq. (2.5) is that it 
respects this equality also for subtracted scattering rates. 

Subtracted rates for processes C, D, E, G, I vanish, and looking at Goldstinos one can easily 
understand why: a single Goldstino diagram contributes, such that no interference terms exist. 
This is not the case for scatterings H and J, where a second Goldstino diagram contributes, 
generated by the quartic Goldstino coupling in eq. (A. 29). (This extra coupling is not present 
for the ghost scatterings in A and B analogous to H and J, as we employ a non-supersymmetric 
gauge without ghostinos). In case of scattering F the subtracted rate vanishes because propor- 
tional tos + t + u = 0. Al/2! factor must be included for the A and F processes that have 
equal initial state particles, and a factor 2 for C, D, G, H that can occur with particles and 
with anti-particles. The total result for the subtracted gravitino production rate is 

where the numerical factor accounts for the difference with respect to the scattering rate com- 
puted in Boltzmann approximation, where 7 ~ aT^/ir^ where a = '^J^^dt |=e/p/167rs^ 
is a constant. The sum runs over the three components U(l)y ® S\J{2)l ® SU(3)c of the 
MSSM gauge group with = {1,2,3}, and Cn = IT^I^ = A^(A^^ - 1) = {0,6,24} and 
C'jy = l^jj P = {11; 21, 48} where J2<s> ^^^^^ over all chiral multiplets. We use the standard 
normalization for hypercharge, where left-handed leptons have Y = —1/2, that differs from the 
SU(5) normalization by a factor a/3/5. All parameters are renormalized at an energy scale 
/i ~ T. 

The next step is computing diagram D: we first need to introduce finite temperature effects. 

3 Finite temperature effects 

We here summarize some well known results from quantum field theory at finite temperature 
that are relevant for our computations: the spectral densities of scalars, fermions and vectors 
that play a role analogous to parton densities in hadron scattering processes. This section 
also contains a few original points: practical formulae for thermal masses that apply to generic 
supersymmetric models, the observation that thermal effects respect supersymmetry a.t E ^ T; 
we explain what qualitatively changes and why we must go beyond the Hard Thermal Loop 
approximation; we discuss a possibly non-standard point of view about the problem of negative 
spectral densities. 
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Figure 6: Dispersion relations at finite temperature in HTL approximation for the components 
within a chiral (left) and vector (right) massless super-multiplet. Thermal effects are supersym- 
metric at k ^ m ^ gT. 



3.1 The Hard Thermal Loop approximation 



Thermal corrections simplify when one restricts the attention to diagrams with ^soff external 
momenta, <C T [6, 5]. This approximation is useful if couplings are small, <C 1, as it 
describes collective phenomena that develop at energies of 0{gT) via simple effective thermal 
Lagrangians. In the rest frame of the plasma, the non-local HTL Lagrangian for scalars (p, 
fermions ip and vectors is [6, 5] 



HTL 



rriglA 



2 I 2 

+ rrip 



ip - m^Tr / 

Jn I 



in K-D' " Jn ""(if ■ ^^'^^ 

where • • • denotes Yukawa or scalar couplings that do not receive HTL corrections; gauge 
couplings receive thermal corrections such that ^rtl is gauge invariant (indeed D denotes the 
usual gauge-covariant derivative); K = (1, fc) is the 'loop' momentum {K^ — ^)'^ In ~ J dfl/in 
denotes angular average. It is performed analytically in the more explicit results in appendices B 
and C 

The key parameters are 'thermal masses' of order m ~ gT. By explicit computation we find 
the following values for thermal masses in an unbroken supersymmetric theory with massless 
chiral $ = (0,0 vector V = (V^, A) superfields:^ 



rrij 



2mi 



T' 



2mi 



9 



(3.2) 



where g is the gauge coupling and A the coupling in the superpotential W = A$$'$". Sum- 
mation over gauge, flavor and any indices is understood. The group factors Cr and are 



^0 is a complex scalar, ^ and A are Weyl fermions. Explicit formulae for thermal masses of bosonic sparticles 
had been given in [7]; we agree with their results. 
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defined as TrT^T^ 



T'^5°'^ (index of tlie representation) and as (T^T^)ij = 0^6^ (quadratic 



Casimir) wliere tlie generators are in tlie representation R. By summing botli over ij and over 
ab one finds that they are related by T^dimG = C^dimi?. Explicit values are = 1/2 and 
Cr = (iV2 _ iy2N for the fundamental of S\J{N) (dimi? = A^, dimG = iV^ - 1), Cy = N for 
the adjoint of SU(A^), and Cr = for a representation of U(l) with charge q. In the MSSM 
with 3 generations and one pair of Higgses one has the following vector thermal masses 
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^V3 - 4 

and the following scalar masses 



rrii 



9y rp2 



ml 



m 



9 



m 



V2 



9lT\ 
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m 



Vi 



9lT^ 



'9 fi'v"' 



9l 



92 + ^ + -X 



rp2 



(3.3) 



(3.4) 



m 



Q 



3^3 ^ ^ 72 4 



ml 



2 Q 2 Q 

5* + 5*- + T 



rp2 



m 



D 



2 2 £/y 
-9I + — 
3^ 18 



rp2 



where the At terms are present only for third generation squarks, and we neglected analogous 
Aft and A,- terms, possibly relevant if tan/? ~ mt/m^,. Squared thermal masses for gauginos, 
higgsinos, quarks and leptons are a factor 2 smaller, as summarized in (3.2). 

We followed the standard convention for thermal masses. Let us recall how they parameter- 
ize thermal dispersion relations ijj{k) where uo and k are the energy and momentum with respect 
to the plasma rest frame. Scalar thermal masses m^ correspond to the relativistic dispersion 
relation uj"^ = k'^ + m'^, see eq. (3.1). For fermions the thermal mass m tells the energy at rest of 
particle and hole (or 'plasmon') excitations, Lj{k = 0) = m, while at large momentum the hole 
disappears'^ and particles have a;^(A; ^ T) ~ fc^ + 2m^. For vectors the thermal mass m tells 
the dispersion relation of transverse polarizations at large momentum, Ci;^(A; ^ T) ~ /c^ + m^, 
while at rest both transverse and longitudinal polarizations have energy uj'^{k = 0) = 2m^/3. 

Therefore, despite the misleading conventional factors 2, eq. (3.2) means that within each 
multiplet, vector or chiral, thermal effects at k ^ T modify in the same way the dispersion 
relation of its bosonic and of its fermionic components. This likely is a consequence of the 
eikonal theorem, that tells that gauge interactions with soft vectors do not depend on the 
particle spin but only on its gauge current (thermal masses physically describe the kinetic 
energy that a particle acquires due to scatterings with the thermal plasma). Fig. 6 shows the 
dispersion relations uj{k) of the particles within chiral and vector multiplets. Particles and 
sparticles have similar dispersion relations, reducing the phase space for gravitino production 
via decays. 



3.2 Full one loop thermal effects 

The HTL approximation holds at momenta and energies k,u <^ T, correctly describing ther- 
mal effects that u ~ gT if g <^ 1 . However, the physically relevant values of gauge 
couplings (especially the strong coupling) are not small enough to justify the use of the HTL 
approximation. We therefore use the full one-loop thermal and quantum corrections to propa- 
gators of scalars, fermions and vectors. Explicit expressions are collected in appendices B and C 

■^More precisely, its residue at the pole is exponentially suppressed by k"^ /m? . The fact that residues Z{k) are 
not constant is one reason why computing the imaginary part of the gravitino propagator in terms of particle 
and sparticle spectral densities is a better formalism than directly computing the gravitino production rate: it 
precisely dictates how all these non- relativistic factors must be taken into account. 
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HTL limit, m ^ T m = T 
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Figure 7: Spectral density of a massless fermion in a thermal plasma, plotted in HTL approxi- 
mation (g ^ 1, i.e. thermal mass m ~ gT <^T) and beyond. Notice the main differences: the 
particle (u > 0) and 'holes' (u < 0) poles develop a finite width, and, more importantly, the 
continuum below the light-cone gets Boltzmann suppressed at k>T. Contours are equispaced 
in log scale. 



(see also [8, 9, 10, 5]) and fig. 7 illustrates (in the case of a massless fermion) the qualitatively 
new effects that arise beyond the HTL limit. 

The most visible effect (although not the most important one) happens at \uj\ > k i.e. 
'above the light cone'. In the g ^ limit particles (and quasi-particles such as fermion 'holes') 
have an infinitesimal width: their dispersion relations are plotted as thin lines in fig. 7a. For 
finite g they get a finite width F (both from T = quantum effects and from thermal effects), 
such that their spectral density gets smeared acquiring the usual bell-like shape. This is why 
a continuum appears also above the light cone in fig. 7b. A well known problem encountered 
in thermal computations is that sometimes thermal effects give F < 0. We therefore included 
the T = contribution, finding that the total F is positive for scalars and fermions. This cure 
does not work for vectors, because the T = contribution to their F can itself be negative 
depending on the gauge choice, see eq. (B.3a). We therefore think that the negativity of F is 
not related to higher-order subtleties in the thermal expansion, but just to gauge invariance. 
It should not affect computations of physical gauge-invariant quantities, provided that one can 
do an exact computation up to some order in the perturbative expansion. The only trouble 
is that in practice it is difficult to achieve this in finite temperature computations. In view of 
this situation, since the would-be poles are anyhow reasonably narrow for the physical values 
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of the coupling that enter our computation, we use for them the HTL approximation^. Notice 
that the HTL approximation correctly describes the position of the poles (i.e. the dispersion 
relations) even at large k>T: poles lie close to the light-cone, |a;| ~ /c, even if ~ 1 [10]. 

The new effect important for our purposes arises at |u;| < k, i.e. 'below the light cone'. 
Quantum effects do not give any contribution to spectral densities here (and more generally 
below the threshold for zero-temperature decays), and the purely thermal contribution is not 
problematic. Even in HTL approximation, thermal effects give non zero spectral densities below 
the light cone: this describes 'Landau damping' i.e. the fact that particles exchange energy with 
the thermal plasma. However the HTL approximation cannot be applied at ~ T (a region 
relevant for us, since g ^ 1), and indeed it misses one key physical fact: at k ^ T spectral 
densities get suppressed by an exponential Boltzmann factor. Indeed the thermally averaged 
coupling of a particle with large momentum A; ^ T is small, since very few of the particles in the 
plasma have the large momentum demanded by energy-momentum conservation. This Boltz- 
mann suppression of the spectral density below the light cone is the main difference between 
fig. 7a (HTL approximation) and fig. 7b (full one loop), and makes the gravitino production 
rate about 50% smaller than what one would find by applying the HTL approximation at all 
momenta, outside its domain of validity p <^T. 



3.3 Vector and gaugino propagators 



We are now going to present how spectral densities are practically used. We need the re- 
summed propagators for the vector with four-momentum K = (/cq, k) and the gaugino with 
four-momentum Q = (go, q) in the loop. We employ non-time ordered propagators (as they 
allow slightly cleaner formulae than imaginary parts of propagators), denoted with a ^ in the 
notation of [5] that we follow. Thermally resummed propagators are denoted with a *: they 
are (see appendices B, C for more details) 



*S<{Q) 



(7o - 7 ■ q)p+{Q) + (7o + 7 ■ q)p-{Q) 



(3.5a) 
(3.5b) 



Some explanations are in order. First, go > or /cq > describes a fermion or a vector in 
the final state, and go < or /co < describes a fermion or a vector in the initial state: this 
convention allows to compactly describe all possible processes. Indeed the factors 



fniko) 
fpiqo) 



1 



e<io/T 



riB if fco > 

-{I + ub) if fco < 

rip if go > 
1 — Up if go < 



(3.6a) 
(3.6b) 



^ This might be not an entirely satisfactory approximation for the pole-pole contribution to gravitino pro- 
duction, as particle and sparticles happen to have similar dispersion relations at fc ~ T, and what matters for 
the phase space is their mass difference. Due to this reason, we will find that the pole-pole contribution is small, 
and it seems unlikely that adding a finite width can change this conclusion. 

Furthermore, one might compensate this approximation by not subtracting modulus squared of s-channel 
diagrams when computing subtracted rates. Since these details have negligible numerical significance, we 
prefect to avoid them. 
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give the usual statistical factors: —n (number of particles in the initial state) or l±n (stimulated 
emission or Pauli-blocking in the final state), where UB^piE) = l/(e'^'/"^=Fl) are the usual Bose- 
Einstein and Fermi-Dirac distributions. 

Second, p+, p_, pr, Pl are the spectral densities for the fermion, fermion pole, transverse 
vectors and longitudinal vectors respectively. As discussed in the previous section, we can keep 
the HTL pole approximation outside the light cone, so that 



p±(Q) = 271 
PlAK) = 27r 



+ pT\Q), (3.7a) 



ZL,T{k) 6{ko - UJL,T{k)) - ZiT^ko + U^L,T{k)) 



+ pTt{Q)- (3.7b) 



In HTL approximation the residues at the poles are given in terms of the pole positions oj±{q) 
and uji^k) as [9, 8, 5] 

^ ujl-q^ ^ ujl{ujl - k^) y ^ ^t{^t - k^) .o on 

^ 2ml ' A;2(A;2 + 2m2,-cui)' 2m^yUjl-{ujl-k^)^' 

These formulae tell that residues for longitudinal and hole excitations are exponentially sup- 
pressed at energies larger than gT: they are low-energy collective phenomena. The continua 
^cont Qj^iy exist below the light cone, at |go| < q and \ko\ < k. The spectral densities satisfy 
sum rules such as 

dqo , r^dko . f-'^^dko 2ml qn 

-P±{Q) = 1, / -7^Pt[K) = 1, / —pl{K) = — — (3.9) 



27r^"^"^ ' ./_oo 27r^^^^ ' .L^ 27r ^ 3P 



and the continuum turns out to contribute ~ (10 -^ 20)% less than the poles. Eq. (3.9) means 
that the number density of longitudinal vectors diverges at A; ^ 0, but this leaves finite gravitino 
rates thanks to the d^k integration factor. In the T = limit uj±{q) = ±q, uji^k) = k and one 
can check that the standard expressions for the propagators are recovered. Notice that pl,t 
have dimensions mass~^, while p± have dimensions mass~^. 

4 Gravitino production rate due to decay effects 

We can now compute the imaginary part of diagram D in fig. 2, and extract from it the gravitino 
production rate. Using the gravitino ^ Goldstino equivalence, eq. (2.4), diagram D is obtained 
from eq. (2.2) by inserting the quadratic parts of the MSSM supercurrent (A. 28) and of its 
divergence (A. 29): 



^(^2) = - E ^^.?^[7^717VA(^) - V2 [(9>.)*(7Va) + (5>.)(7'^7^ei^)] 



N=l 

3 



N=l 

where runs over the three factors of the MSSM gauge group, and Fu^^ here stands for the 
linearized part of the corresponding field strength. We ignored soft-breaking squared masses of 
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scalars, as they have higher dimension than gaugino masses M^v- The contribution to 11'^ from 
diagram D is 



n< 



1 



2 



Tr((9-5)(2)f(9-5)(2))T 



(4.1) 



Sm'i 



'3/2 



3 



32M2 



1 




) 



(4.2) 



We now explain how eq. (4.2) is obtained. The Goldstino part, proportional to M^/Smg^g; 
is straightforward. We emphasize that the divergence of the supercurrent is evaluated before 
evaluating its thermal matrix element. Indeed, while thermal masses naively look like SUSY- 
breaking terms of order g'^T, they actually do not contribute to d^S^, and a mistake about this 
issue would make the Goldstino rate qualitatively wrong [11, 12, 13] (see also appendix A). 
Indeed, despite the nice formalism employed to compute them (periodic and anti-periodic 
boundary conditions in imaginary time for bosons and fermions respectively), thermal effects 
just are one particular background: no background affects the operator equations of motion, 
such that a supercurrent which is conserved at T = remains conserved at finite T.^ 

For the remaining massless gravitino part, we insert the explicit value of the gravitino 
polarization tensor (2.5) and get two terms of the form 



where ^^(2) = T/^'S'^^^. We can now perform simplifications that only employ the known Dirac- 
matrix structure of 5*1^^: 

• The vector /gaugino contributions obey ,^(2) = (thanks to 7/x[7a, 7/3]7'^ = 0), such that 
only the first term of eq. (4.3) contributes. It is reduced to the same operator On using 
the '~ifj,'~ial^ = —27a identity and taking into account that the thermally corrected gluino 
propagator has the same 7-matrix structure as the massless propagator. This leads to 
the 1 + M|r/3m3^2 prefactor in eq. (4.2). 

• The quark/squark contributions vanish, thanks to a cancellation between the two terms 
in eq. (4.3). Indeed, by applying the ■j^'jal'^ = —'^'Ja identity (one time in the first 
term, and two times in the second term) both terms reduce to the matrix element 
Tr{^R{dip)f{dip*)C,L), with opposite coefficients. 

We don't know if there is some deeper reason dictating these cancellations such that the 
full result is controlled by the thermal matrix element of the operator times the prefactor 
1 + Mjf /Sm^^^. A general proof of this result would allow to get the full production rate from 
the simple Goltstino rate according to eq. (4.2). 

For completeness we mention that we have studied thermal corrections to the On operators 
in HTL approximation. As well known gauge vertices g receive very large thermal corrections 

^Although this is not relevant for us, we can be more precise: thermal effects spontaneously break super- 
symmetry in the visible sector, and the associated thermal Goldstino mode was identified with a particular 
collective excitation [14]. The conservation of the supercurrent at finite temperature is therefore analogous to 
how electrowcak gauge currents remain conserved despite the Higgs vev. However, since the thermal Goldstino 
is a low energy phenomenon, we don't know how to extend it to write an explicit conserved supercurrent that 
also holds at energies E ^ T. 



TY(^f,/5f2) )t + ^Tr(^(2)f ^(2) )t 



(4.3) 
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of order g{l + g'^T'^/k'^) where k <^T (HTL approximation) is some external momentum: their 
presence would be problematic, as they seem to describe infra-red divergent effects (see e.g. 
section 10.3 of [5]). In the case of gauge vertices these corrections are demanded by gauge in- 
variance: different diagrams combine such that =^htl of eq. (3.1) contains the gauge-covariant 
derivative D. On the contrary Yukawa couplings do not receive these problematic HTL cor- 
rections. We verified that the Goldstino vertex x^n does not receive any HTL correction.^ 
Beyond the HTL limit there will be corrections suppressed by powers of g/ir, that we can ignore. 



4.1 Gravitino propagator 

We now restart from eq. (4.2) and explicitly compute the imaginary part of the gravitino 
propagator with four-momentum P = [po, p) = K + Q, summed over its polarizations: 

(4.4) 

where N = {1,2,3} runs over the three factors of the SM gauge group with un = {1,3,8} 
vectors; Mj^ are the gaugino masses at zero temperature (renormalized at some scale around 
T). Inserting the explicit parameterization 

^^=(^0,^,0,0), Q = {qo,q cos 6q,q sin 6q), P = {p,p cos Bp, p sin 6p) 

for the vector, gaugino and gravitino four-momenta respectively one finds 

r2 \ r 



(27r) 



fB{ko)fF{qo) X 



;V=i V "^""3/2/ ^"Pl 

PL{K)p.{Q)e cos' + pL{K)p^{Q)e sin' + 



+PTiK)p+{Q)^{k' + kQ){l + COS 6p COS 6g) - 2kko{cos 6p + COS 6g) ) + (4.5) 
+Pt{K)p_{Q) I {k' + /co)(l — cos9pCos6g) — 2kkQ{cos 6p — cos 9q) 



To compute the total rate 'Jd using eq. (2.1) it is convenient to multiply by 1 = / d^Q S^{K 
P — Q),'^ perform the non-trivial angular integrations over 6p and 6q, obtaining 

Id 



2(27r)3M|i 



N=l ^ "^"^3/2 / 



^The basic reason is the following. Since the Goldstino vertex has dimension 5, by dimensional analysis it 
receives gauge corrections of order J d'^K{K^)^/[K'^{K + PiY{K + ^2)^] where {K^)fj. denotes any vector 
formed with 3 powers of K: it necessarily contains the combination A'^, that, as explained in [5], docs not lead 
to HTL vertices. 

^This step also allows to see that the seemingly esoteric expression is actually equivalent to what one would 
naively guess from the kinetic theory, if spectral densities are treated like parton densities 

377^ = E / ^^^,P±{Q)PL.T{K)\s^?{'2^f5\P ±Q±K)- (statistical factors) 
dV dt dP ^ Jqo,ko>o (27r)^ [2n)-' 

where the sum is over all polarizations, gauge indices, gravitino production processes with amplitudes =2/. As 
discussed around eq. (3.6), the factors /^(fco) a-nd /f(9o) reproduce the usual statistical factors, 1 ± n or —n. 
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where 



N 



T 



-6 



dqo dko kfB{h)fF{qo) x 



PL{K)p.{Q){p - qY[{p + qY - k'] + pL{K)p^{Q){p + qY[e - (p - qY] + 



+Pt{K)p^{Q){{p + qY - k') (1 + k'Jk'Yk' + {p- qY) - Ak,{p - g) + (4.7) 



+PT{K)p_{Q){k' -{p- qY) (1 + k'je){k' + {p + qY) - Akoip + q) 



The dimensionless coefficients are positive: each term in the square brackets is positive 
in the allowed region, except pr that becomes positive after being multiplied by fsiko). The 
integration range is restricted by momentum conservation, p + k + q = 0, i.e. \k — q\ < p = 
ko + qo < k + q: any side of a triangle cannot be longer than the sum of the other two or shorter 
than their difference. 

The last two equations generalize eq. (38) of [2], who considered the vector /gaugino loop 
in the limit of hard gravitino and soft vector (small ko <^ P,T, such that fsiko) — T/k^) and 
neglected the gaugino thermal mass (i.e. p_ ~ and p+ ~ 2n6{qo — q) such that q = qo = p—ko). 



4.2 Decay contribution to the gravitino production rate 

In conclusion, the decay contribution to the gravitino production rate per space-time volume is 
given by eq. (4.6). The coefficients have to be evaluated numerically. We approximate the 
spectral densities outside the light cone as 5-function poles (using full expressions they would be 
narrow bells, making numerical integration difficult for our limited computing power), such that 
we have four types of contributions: pole-pole, continuum-continuum, (vector pole)- (gaugino 
continuum) and (vector continuum)- (gaugino pole). A vector can be either longitudinal or 
transverse, in the initial state or in the final state, and similarly for the gaugino. 

The resulting coefficients /at depend on the gauge couplings and on the content of matter 
charged under the given gauge group; in the MSSM it is convenient to parametrize them as 
functions of the thermal vector masses my^ listed in eq. (3.3): 

fN ^ InC^). (4.8) 

For example ms ~ 1.3T for the gluon at T ^ 10^ GeV. The functions /at are plotted in fig. 1. 
In HTL approximation there would be a unique iV-independent function /, and the functions 
/tv turn out to be somewhat different, depending on the relative amount of vector and chiral 
multiplets present within each group. In non-minimal models with more chiral multiplets than 
in the MSSM, one would have to add their extra contributions to vector thermal masses, and 
to slightly revise the functions /at. 

Finally, let us try to discuss the accuracy of our result. Thermal corrections to the pressure 
have been computed up to high orders in [15]: these computations can be used to see how 
convergent the perturbative expansion is in practice. In the favorable limit Nj 3> Nc (where 
Nf is the number of flavors and N^. is the number of colors) the perturbative expansion for the 
pressure remains accurate up to my /T 1 (fig. 1 of [15]). This presumably also applies to our 
case, as SUSY-QCD has a set of fermions and scalars that give the same contribution to the 
gluon thermal mass as Nf = 21 flavors in the fundamental. 
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Figure 8: Top scatterings. Feynman diagrams contributing to QU ^ . The • denotes 
the top Yukawa coupling. 




Figure 9: Top scatterings. Feynman diagrams contributing to QU ^ p.H . The • denotes a 
coupling proportional to Xt, and o its A-term. 



Furthermore, AdS/CFT techniques should allow to compute the large coupling limit of the 
gravitino emission rate in some supersymmetric theory, maybe not unrealistically different from 
SUSY-QCD. This could be done analogously to how [16] used AdS/CFT to compute the photon 
emission rate from strongly coupled Af = 4 SYM in the large Nc limit. By analogy, we expect 
that at strong coupling the gravitino rate functions will have a finite limit, A^c-independent 
up to 1/A"c corrections. 

5 Production of gravitinos due to the top Yukawa 

Previous works considered gravitino production due to the g^, g2 and gy gauge couplings; 
the top quark Yukawa, QUH, also has a sizable coupling A^. There are two main kind of 
scattering processes: 

a) Scatterings involving fermions only, such as QU — > "^H: fig. 8 shows the relevant Feyn- 
man diagrams. They would populate only the spin 3/2 component of the gravitino, as 
only dimension-2 soft terms enter these diagrams, so that Goldstinos are not produced. 
However, an explicit computation shows that the dominant contribution, of order T^/Mpj, 
vanishes. 

b) Scatterings involving two fermions and two scalars, such as QU — > "^H: fig. 9 shows the 
relevant Feynman diagrams. The first diagram involves At, the dimension-1 A-term of 
the top Yukawa coupling, and populates the spin 1/2 component of the gravitino. The 
other three diagrams populate the spin 3/2 component. 
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Figure 10: Evolution of the gravitino ahundancy nj s in units of '-f/H^s at T 
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The total result is: 



|i2/(top scatterings) 

all 



72 



(5.1) 



where s 
rate is 



(Pi + P2)^ is the usual kinematical variable. The corresponding gravitino production 



Itop 



1.30 



(5.2) 



where the numerical factor 1.30 is the correction due to the Fermi-Dirac and Bose-Einstein 
factors with respect to the Boltzmann approximation. 

In the language of previous sections, eq. (5.2) is the scattering contribution. We now ex- 
plain why it also is our total result. First, it happens to be infra-red convergent: the potentially 
divergent contributions given by the modulus squared of the t-channel and ^-channel diagrams 
in fig. 9 actually vanish. Therefore, unlike in the case of gauge scatterings, the inclusion of 
thermal masses is not necessary for obtaining a finite result. Furthermore, including thermal 
effects along the lines of the previous sections does not affect the final result. Indeed the top 
Yukawa coupling gives a thermal mass for top, stops (and higgs and higgsinos): the resulting 
quark/squark/gravitino (and higgs/higgsino/gravitino) decay rates have been computed in sec- 
tion 4 for generic thermal masses, and vanish. Consistency requires that the subtracted top 
scattering rate equals the total scattering rate of eq. (5.2), and indeed the subtracted terms are 
the modulus squared of the t-channel and w-channel diagrams in fig. 9, which vanish. 

Again, all these cancellations have a simple interpretation: they are the ones needed such 
that the production rate for the spin 3/2 components of the gravitino equals the production 
rate for the spin 1/2 Goldstino components, up to the prefactor in eq. (5.1). Indeed, using the 
gravitino/Goldstino equivalence, the Goldstino production rate can be equivalently computed 
from one single diagram that only involves the single quartic Goldstino coupling, 

AfXt Goldstino (higgsino squark squark* + quark squark* higgs + h.c), 

such that decay contributions and subtracted scattering rates simply do not exist for the Gold- 
stino. 
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6 Boltzmann equations with reheating 



We here compute the gravitino abundance by integrating the relevant Boltzmann equations. 
While previous works ignored the history of the universe prior to its reheating (from the point 
of view of computing gravitino production this in practice amounts to assuming that the Big 
Bang started at the maximal temperature Trh), we here follow the standard definition of the 
reheating temperature Trh, where MSSM particles are progressively reheated by the energy 
released by some non-relativistic energy density p^, which could describe e.g. an oscillating 
inflaton field, or some non-relativistic particle decaying into MSSM particles.® In both cases 
the relevant Boltzmann equations are 



P0 + 3Hp^ = 
PR + AH PR = r 

723/2 + 3if7i3/2 = 



7 



(6.1) 



where 713/2 is the gravitino number density summed over its polarizations, a dot denotes d/dt, 
H = R/ R = ^J^Tt{p^ + pr)/3/Mpi is the expansion factor, pr = n'^g^T'^/SO is the energy 
density of MSSM radiation at temperature T (with g^, = 228.75, up to Oi^g"^) corrections, 
and up to adding right-handed neutrinos), and parameterizes the decay width of p^. The 
reheating temperature Trh is defined in terms of as the temperature at which [18] 



I.e. 



RH 



45 



A-K^g^ 



.1/4 



(6.2) 



It is convenient to rewrite equations (6.1) in terms of Y{z), where z = T^ji/T and Y = n^/2/s, 
with s = 4:Pr/3T being the MSSM entropy density. Following [19] one gets 



HZz 



dz 
dY 



—3Hp^ — T^p^ 



(6.3) 



sHZz— = 3sH(Z -1)Y + -f 
dz 



where 



PR 



F 



1 _ ^i^. (6.4) 
AprH AHpR ^ > 

To clarify the physical meaning of Trh, we emphasize that Trh is not the maximal temperature; 
however what happens at T ^ Trh gets diluted by the entropy release described by the Z — 1 
factor (Z~3/8atT:» Trh and Z ~ 1 at T -C Trh)- In our case 7(T) oc and the solution 
is 



Y{T < Ti 



RHj 



7 

Hs 



0.745 



7 



T=Tr 



HrS 



6.11 10 



-12 



RH 



7|t=Trh 



T=Trh 



1010 GeV T^jj/M] 



(6.5) 



Fig. 12 shows our results for the dimensionless order-one combination 7/(T^/ Mpj) that appears 
in (6.5). Notice that the gravitino abundance is proportional to it, and to Trh^ the large power 
7(T) oc gets almost compensated by cosmological factors. 



^ Alternatively, some of the flat directions present in the MSSM supersymmetric potential might develop 
large vevs during inflation. There is a debate in the literature whether such condensates can be sufficiently 
long-lived to affect reheating [f7]. For simplicity, we here do not consider these possible but model-depenent 
phenomena. 
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In previous analyses p<j(, was ignored and the 'instantaneous reheating' approximation was 
used, which amounts to start the Big Bang from a maximal temperature T = T^u'- this gives 
a slightly larger gravitino abundance 



Y{T « Trh) ^ 



(6.6) 

T=Trh 



Fig. 10 illustrates the different evolution of Y (in arbitrary units) between the two cases. 

Within the standard ACDM cosmological model, present data demand a DM energy density 
^Buh"^ = 0.110±0.006; if DM are non relativistic particles with mass M ^ keV this corresponds 
to Fdm = (0.40 ± 0.02) eV/M [20]. One has Y = Fdm if gravitinos are the observed DM. 
Equivalently, one can compute the present gravitino mass density in terms of their relative 
entropy Y as 

^'^'^ - -JjW - ^ - °-°°^^^GeV10iOGeVT^^/M|/ ^^'^^ 

where pcr = SHqMpJStc is the critical energy density, Hq = 100 /ikm/sec-Mpc is the Hubble 
constant, the present entropy density is Sq = 2'7i'^g^sTQ /4:5 with g^g = 43/11 and Tq = 2.725 K. 
Fig. 11 compares the regions where the thermal gravitino abundance equals the DM abundance 
with the regions compatible with standard thermal leptogenesis [24], as computed in [19] with 
the same definition of the reheating process. This plot ignores all model dependent issues, 
including who is the LSP and the NLSP. Let us briefly summarize these issues [1, 21, 22, 23]. 

• If the gravitino is the stable LSP then 

— if m3/2 ^ keV the gravitino behaves as cold dark matter and its energy density can 
be at most equal to the DM density. 

— a somewhat stronger constraint applies if the gravitino is lighter, Tq -C 777,3/2 < keV, 
and consequently behaves as warm dark matter or radiation. The Goldstino compo- 
nent of such a light gravitino thermalizes (unless Trh is as low as possible), so that 
this scenario is severely constrained: assuming that the Goldstino was in thermal 
equilibrium when g^, ~ 100, present data demand 7773/2 ^ 16 eV [22]. 

An additional contribution to the gravitino energy density, ^^3/2'^'^ — fiNLSp7?73/2/7?7NLSP 
(having neglected entropy production, which typically is an excellent approximation) is 
generated by NLSP decays, with mass 777nlsp and mass density I^nlsp after their freeze- 
out. Weak scale sparticles give ^nlsp ~ 1? such that this extra contribution is significant 
if 7773/2 is not much smaller than m^LSP- 

• If heavier than the LSP (which possibly has mass tt^lsp ~ 100 GeV), the gravitino gravi- 
tationally decays into the LSP and some SM particles: 

— If 7773/2 ^ lOTeV the gravitino decays before BBN, generating a contribution to 
the LSP energy density, flhsp — 7?7lsp ^3/2/7^73/2 [1] (we neglected the entropy in 
gravitinos). 

— A lighter gravitino decays during or after BBN, damaging nucleosynthesis. The 
resulting bound on ^3/2 depends on which SM particles are produced by gravitino 
decays, and typically is some orders of magnitude stronger than the DM bound 
^3/2^^^0.1 [21]. Very late gravitino decay into photons would also distort the 
CMB energy spectrum. 
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Gravitino mass in GeV 



Figure 11: The bands show the region 
where the thermal gravitino abundance 
equals the DM abundance (3a regions), as- 
suming unified gaugino masses with mi/2 = 
150 GeV (roughly the minimal value al- 
lowed by present data) or mi/2 = 1 TeV 
at the unification scale, and negligible At. 
Model- dependent issues are here ignored, 
including who is the LSP and the NLSP. 
Successful thermal leptogenesis with zero 
initial right-handed neutrino abundance is 
not possible within the gray band [24, 19]. 



We recall that we computed thermal production of gravitinos from MSSM particles at tempera- 
tures Trh ^ "^3/2, "^soft- The true physics might be different. For example, the messenger fields 
with mass Mgm employed by gauge mediation models might be so light that they get thermal- 
ized (together with the hidden sector) and later decay back to MSSM particles, leaving a ther- 
malized Goldstino. Eq. (6.3) shows that this phenomenon is dominant if Mqm ^ (10-i- 100)Trh, 
as the gravitino abundance gets washed-out as y oc during reheating at T> Trh- 



7 Conclusions 

Previous computations of the thermal gravitino production rate [2, 3] were performed at lead- 
ing order in small gauge couplings, finding a rate of the form 7 oc g'^lnl/g, which behaves 
unphysically when extrapolated to the true MSSM values of the gauge couplings, g ^ I (see 
fig. 1). We improved on these results in the following ways: 

1. We included gravitino production via gluon — gluino + gravitino and other decays: these 
effects first arise at higher order in g (the phase space is opened by thermal masses), but 
are enhanced with respect to scattering processes by a phase-space vr^ factor, typical of 
3-body vs 4-body rates. The gravitino production rate becomes about twice larger, or 
more if M3 > Mi_2 > "^3/2- 

2. We added production processes induced by the top quark Yukawa coupling. This enhances 
the gravitino production rate by almost 10% or more if At is bigger than gaugino masses. 

3. Finally, we computed the gravitino abundance replacing the instant reheating approxima- 
tion with the standard definition of the reheating process, where Trh is not the maximal 
temperature but defines the temperature at which infiaton decay ends, ceasing to re- 
lease entropy. This improvement decreases the gravitino abundance by 25% and allows a 
precise comparison with leptogenesis [24], where reheating was included in [19]. 

Our result for the gravitino production rate is 

7 = 7D + 7s"' + 7top. (7.1) 
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Figure 12: Production rate 7 = dN/dV dt in units ofT^/Mp^ for the spin-3/2 gravitino compo- 
nents in the MSSM. The upper curve is the total rate, and the other continuous curves show the 
contributions from {qs, g2, Qy, ^t} interactions (summed over decay and scattering processes). 
The production rate for the Goldstino spin- 1/2 components is obtained by multiplying these 
four contributions times {M|, M|, Mf, Af}/3m^^2 respectively. The dotted curve show previous 
results from [2] and [3]. 

where the decay rate 7d (which dominates the total rate) is given in eq. (4.6), the subtracted 
scattering rate 7™^ in eq. (2.6), and the rate induced by the top Yukawa couphng in eq. (5.2). 
Fig. 12 summarizes our results, showing the value of the dimensionless combination 7/(T^/Mpj) 
(as well as the values of the single gauge and top contributions to it) which determines the 
gravitino abundance as in (6.5). In this plot we assumed = 0.7 and a unified a = 1/24, 
renormalized at the scale Mgut = 2 ■ 10^^ GeV. 

Accessory results scattered through the paper include: a clean precise re-derivation of grav- 
itino couplings; expressions for thermal masses in a generic supersymmetric theory; the observa- 
tion that they respect supersymmetry at energy much larger than the temperature; a collection 
of formula for thermal corrections to vectors (including correct imaginary parts) beyond the 
Hard Thermal Loop (HTL) approximation; a possibly non-standard discussion of the physical 
meaning of negative spectral densities; a technique that allows to deal with Coulomb-like infra- 
red divergences without introducing an arbitrary splitting scale /c* that satisfies the problematic 
condition gT -C fc* <C T. 

A curious result simplified our computation: the differential production rates for the spin 
3/2 and for the spin 1/2 (Goldstino) gravitino component are equal up to a universal prefactor, 
despite that Goldstino couphngs (to dimension-1 SUSY-breaking soft terms in the supercurrent) 
apparently are much simpler than the spin-3/2 couplings (gravitational, to the supersymmetric 
supercurrent). This equality holds thanks to various cancellations, such that various troubling 
contributions automatically drop out from our computation. 

Acknowledgements We thank R. Barbieri, C. Scrucca, G. Giudice, A. Notari, G. Moore. A.S. 
thanks R. Rattazzi who suggested that we look into the higher-order corrections to the gravitino 
production rate because this problem involves conceptual issues of academic interest. On the contrary 
the 100% enhancement with respect to previous results is relevant for phenomenology, and we could 
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compute it without understanding why we don't need to understand how to write the supercurrent at 
finite temperature. 



A Gravitino propagator and couplings 

We here derive the needed gravitino propagator and couphngs, both generic and speciaUzed for the 
MSSM. The results contained in this section are not new, but they are useful for two reasons. First, 
because all factors (z, 75, P^, etc) and subtleties should be right, and are relevant for satisfying the con- 
sistency checks that we performed in our subsequent computations. Second, because we recomputed 
relevant gravitino properties in a way that we consider simpler than in previous literature: we proceed 
directly, without using Noether and supergravity techniques, which are unnecessarily cumbersome 
for our purposes. We use the standard Weyl spinor and 7-matrix conventions corresponding to the 

signature (H ), see e.g. [25]. The phase of the gauginos is chosen such that gaugino couplings to 

matter are real. We assume a Minkowski background i.e. we neglect the small cosmological constant. 

A.l The gravitino Lagrangian 

The gravitino is the gauge field associated with local supersymmetry, and becomes massive by means of 
a super-Higgs mechanism: 'eating' the massless Goldstino fermion arising when global supersymmetry 
is spontaneously broken. This is analogous to a gauge vector that becomes massive via the usual Higgs 
mechanism, so that we start recalling some general properties of this well known simpler case, and 
this analogy will later allow us to derive gravitino properties following the same logic. 

Paradigmatic digression 

We thus consider a U(l) gauge symmetry broken by the vev of a charged scalar field H. In the limit of 
vanishing gauge coupling, a massless Goldstone x appears in the expansion of H around the minimum, 
H = v + ix, and x transforms (under a U(l) rotation with infinitesimal angle e) as 6x = v£. The total 
U(l) current is given by 

J^. = Jr'-^d^X, (A.l) 

where J™'^* is the U(l) current of the rest of the theory (e.g. fermionic matter). The total U(l) current 
is conserved: 

This is the case if and only if the Lagrangian contains the coupling 

^D^X-d.jp'. (A.3) 

This is sometimes known as Goldberger-Treiman relation, and shows that Goldstone interactions are 
predicted in terms of nonconservation of the matter symmetry current induced in the process of 
symmetry breaking. 

When the U(l) symmetry is gauged, the total gauge-invariant Lagrangian is 

--/lu + \{d,X - vA,f + {A, - \d,x)Jr' + ^mat. (A.4) 

We can fix the unitary gauge by setting x to zero or, equivalently, by redefining A'^ = — d^xh- 
The second term in (A.4) becomes a mass term for A!^. Notice that while A^ couples to the total 
conserved current given by eq. (A.l), the massive vector A!^ couples to J™^*. 
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The last thing that we want to recall concerns production of massive gauge bosons at high energy 
E ^ V. The effective Lagrangian appropriate for this situation is obtained from (A. 4) by keeping 
terms with the highest number of derivatives in and Xi ^-nd is given by 

^HE = -\fI, + \{d,x? + A.J^-'^' + \x ■ 5^ + =^mat. (A.5) 

We thus see that the total cross section can be approximated (up to terms suppressed by v/E) by the 
sum of massless gauge boson production plus production of Goldstones with coupling (A. 3): 

a{A'^) = a{A^)+a{x). (A.6) 

This statement is called the equivalence theorem. It can also be deduced (in a less transparent 
way) from the fact that the physical state projector for the gauge boson of mass m takes the form 
—Qiiv + k^k^/rn^. Notice that, while to get the correct Goldstone production rate it is crucial to 
take current nonconservation into account, in computing cr{A^) we can actually assume that J™^* is 
conserved. 

The main points of the above discussion — the form of the total current, the Goldberger-Treiman 
relation, the fact that the massive gauge boson couples to the same matter current, and the equivalence 
theorem — will find their analogues in the gravitino case. 



Goldstino interaction 

We now repeat the steps in the previous section in the case of supersymmetry, under which the 
Goldstino x transforms as 6x = V^Fe, where -F is a supersymmetry-breaking vev^ and e is the 
supersymmetric parameter. The supercurrent is 

S'^ = S>:,^ + zV2F^^X (A.7) 

where the apex 'vis' signals that we are interested in theories consisting of a visible and a hidden 
sector. Supersymmetry is broken spontaneously in a heavy hidden sector, and its low energy remnant 
is the Goldstino field: all other hidden sector fields can be integrated out, if one is interested in energies 
below the messenger scale. The full supercurrent is conserved: 

d^s'' = a^s^i, + iV2Fd X = 0. (a.s) 

The vanishing of (A. 8) gives the equation of motion for X) and consequently implies the following 
Goldstino Lagrangian: 

■5<;oidstmo = ^xi&X- "^71^^ ^f^^vis (^-9) 

where • • • indicates couplings involving two or more Goldstinos, not needed in our computation. 



Massless gravitino 

In the supersymmetric limit, the massless gravitino is described by a Majorana Rarita-Schwinger field 
-0^ with Lagrangian 

^ = Ae'^'^P-i^^j^ludp^l^. (A.IO) 

invariant under the gauge SUSY transformations with parameter e: Sip^ = —2M-p\d^e. Here Mpi = 
Mpi/Svr = 2.4 10^^ GeV is the reduced Planck mass. The variation of the matter action defines the 
Majorana supercurrent as 

<55^atter = / d^X S^'d^E. (A.ll) 



^Although we choose the notations and normalizations which arc standard for _F-tcrm supersymmetry break- 
ing, the result apply to any combination of F and D-tcrm breaking. 
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Demanding that the fuh action is invariant to zeroth order in M-p^ one obtains how the massless 
gravitino interacts with the supercurrent 



^int = -^i^nS^". (A.12) 
2Mpi ^ ^ ' 



Super-Higgs mechanism 



We win now follow how the massless gravitino eats the Goldstino, getting a mass via the super-Higgs 
mechanism. First of all, the gauge-invariant action for the goldstino-gravitino system is [26, 27]: 



(A.13) 



It contains a gravitino-goldstino mixing mass term, that agrees with the form of the supercurrent, 
eq. (A. 7). Indeed, this Lagrangian is invariant under the local field transformations: 

Sij^ = -Mpi{2d^e + imy2l^e), 
5x = V2Fe, 

provided that the gravitino mass the SUSY breaking vev F are related as 

F 

\/3A?pi ■ 



my2 = ^^7^- (A.14) 



(The derivation above used the flat space assumption). Introducing the gravitino mass required a 
deformation of the supersymmetric transformation of the gravitino, and the gravitino interaction 
term with matter (A.12) is no longer invariant. To restore invariance, we must add to the Lagrangian 
the term 

^3/2 11 



Xi$ = -xi$ $ = -fas'"- (A.15) 

2V2F^ 2MpiV6 ™ ^ ^ 

It may seem surprising at first to find this new coupling of Goldstino to the supercurrent in addition 
to the one in (A. 9). However, there is no contradiction since the new term vanishes as gravity is 
decoupled. 

We can now choose the unitary gauge x = or equivalently define^'' 

1 P^I'X 



such that the whole Lagrangian describing gravitino, goldstino, and their interaction with matter can 
be rewritten as [27] 

This is the Lagrangian describing the massive gravitino We see that it couples to S^-^^. 
Equivalence theorem 

The Lagrangian (A. 17) could be used to study production of massive gravitinos at any energy below 
the messenger scale. In this paper, we are interested in energies much bigger than 7113/2 
sparticle masses. A simpler effective Lagrangian appropriate for these energies can be derived by 



10 



This notation reflects the fact that ^'^ is 'bigger' than V'/i since it contains more degrees of freedom. 
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noticing that the mass terms and mixings between and x in (A. 13) can be neglected. Thus we can 
study production of massless gravitinos and Goldstinos x coupled to the visible sector by 

This is the analogue of the previously mentioned equivalence theorem for production of gauge bosons 
at energies much larger than their masses. In the massive gauge boson case, the equivalence theorem 
could also be derived from the form of the physical state projector of the massive gauge boson. Below 
we will see that an analogous derivation can also be given for the massive gravitino case. Just as in 
the gauge boson case, we can assume that the supercurrent in the coupling S'^-^^'tp^ is conserved; in the 
Goldstino coupling the current nonconservation is of course crucial and has to be taken into account. 

A further simplification concerns the Goldstino coupling (A. 15): as we explain below this coupling 
is irrelevant in MSSM at energies much above the //-term due to approximate scale- invariance. Thus 

It is instructive to compare the relative importance of the two terms in (A. 19) for the total 
production rate. Since the divergence of the supercurrent will be proportional to the soft-breaking 
masses (see below), the effective coupling in the second term is rrisoft/F ~ 1/Mmess- Thus the two 
terms are equally important if SUSY breaking is mediated by gravity, Mmess ~ Afpi- If instead 
-^mess ^ -^Pb like in gauge-mediation models, the Goldstino term dominates [1]. 



A. 2 The gravitino propagator and polarization tensor 

Massless gravitino 

Since the Lagrangian is invariant under local supersymmetry, the same physics can be described by 
different choices of gravitino propagators and polarization tensors. Analogously to the vector case, 
the sum over the two physical transverse polarizations is (see [28]) 

^ E = - ^^M7.f ^ - + P,P.Q (A.20) 

where U is an arbitrary 4-velocity that defines a preferred reference frame used to define what 'trans- 
verse' means, and = [2C/'^(C/ • P) — P^\/2(U ■ P)'^. Gauge-invariant observables do not depend on 
the choice of U. 

As usual, local gauge invariance allows to define more convenient gauge choices. For example, 
one can impose the gauge-fixing condition F = 7^^^ = (in this gauge one also has d^^^ = as a 
consequence of the equations of motion). To derive the propagator, it is best to consider an analogue 
of the ^-gauge by adding to the Lagrangian the gauge-fixing term Fid F/^. Then the kinetic operator 
is invertible, and the gravitino propagator is [29] 

with n^/2 = i7W-(2 + 0^^. (A.21) 



3 /2 

As usual, H/Ju is also the projector to be used when the massless gravitino production rate is summed 
over the gravitino polarizations. The dependence on the gauge-fixing parameter ^ is irrelevant because 
the massless gravitino couples to the conserved supercurrent, Pfj,S^ = 0. For the simplest choice = —2 
one has 

The last two terms do not contribute, again because the supercurrent is conserved. For our later 
computation we will choose 

ii%' = -\i^.riu-rvt.u. (A.23) 



= 7/. = - o7/.f 7. - rv^.. + 7^.P. + P^.1.- (A.22) 
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Massive gravitino and the equivalence theorem 

The massive gravitino is described by the Lagrangian (A. 17). The mass term breaks the gauge 
symmetry present in the massless case. The equations of motion coming from the free part of (A. 17) 
imply 

7^^'^ = 0, d^^f, = 0, (i^ - 7713/2)^-^ = 0. (A.24) 

In the massless case the first two equations could have been imposed as gauge-fixing conditions. The 
resulting propagator is [28] n^,y/(P^ — rriy^ + ie) where 

n^. = -{f + ms/2) \g^.u - ^ ) + if - U + . (A.25) 

V "^3/2/ "13/2/ V 

Again 11^,^ is also the polarization tensor to be used when summing over all physical polarizations: 
n^j. = J2i=±l,±'^ One can check that (A.25) is consistent with (A.24). 

In this paper we are interested in production of ultrarelativistic gravitinos. To study this limit, 
we expand (A.25) in powers of 

+ (-9^.u^ - l7^.flu + ^rP-" + ^t"^'^) + " 5^") ^3/2- (A.26) 

If the supercurrent to which the gravitino couples is conserved, the terms singular in give no 

contribution. The last term vanishes for ~^ 0- The term that does not depend on 7713/2 differs 
from the massless gravitino projector (A. 22) by 7^f^7iy/6, modulo irrelevant terms proportional to 

or Py. 

Thus in the limit 70,3/2 ^ we not only recover the massless gravitino, but also get an additional 
massless spin 1/2 fermion which couples to l/-v/6 times the 'trace' of the supercurrent, $ = ^^xS^- 
This is akin to the van Dam-Veltman-Zakharov discontinuity [30] encountered when adding to the 
graviton a Fierz-Pauli mass term 777^: the limit 77ig then describes the usual massless graviton 
plus a scalar coupled to the trace of the energy-momentum tensor T'^. In our case this 'discontinuity' 
is entirely expected and is consistent with the equivalence theorem as expressed by eq. (A. 18): the 
extra spin 1/2 fermion is nothing but the Goldstino. 

When we take soft SUSY-breaking into account, the supercurrent is no longer conserved. The first 
term in (A.26) can then be interpreted as corresponding to the Goldstino production due to the last 
term in (A. 18) (the coefficient agrees as one checks using (A. 14)). In this derivation of the equivalence 
theorem it is non-obvious that the terms in (A.25) proportional to 'ni~J^ should cancel, as is required 
for full agreement with (A. 18). However this cancellation does happen, as verified in the explicit 
computations needed for this paper. 

A. 3 MSSM supercurrent at zero temperature 
Gravitino couplings 

In a generic renormalizable SUSY gauge theory with vector supermultiplets (^^, A") and matter chiral 
supermultiplets = ((/>«, ^i) and superpotential W the Weyl part s'^ of the Majorana supercurrent 
5^ = [s^^s^) is (see e.g. [31] or explicitly compute it) 

s^^ = -V2 [{D'^<p,na,a^^,) + iW*{cf>*)a'^^,] - ^F^ia'^a^a'^X'^) - ig{<l>m<Pj){a^^\-) (A.27) 

where Df- = 6ijda+igA'^"'T^- is the gauge-covariant derivative. The first two terms are the supercurrent 
of the Wess-Zumino model and of the SUSY gauge theory without matter. The third term is a 
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correction which arises as a result of couphng between the two. (With the Noether formahsm it would 
arise because the Lagrangian is supersymmetric up to a total derivative). In 4-component notation it 
becomes^^: 

S>^ = -V2 [(i?>.)*(7'^7^a) + P>0(7"7^Cy - iWMl^a " iH^*('/'*)7^Cy (A.28) 

where we introduced Majorana spinors and (A", A")"^ which by abuse of notation we denoted 

again and A". As usual = Pl^,, S,R = PrCi = dW/dcpi and the index i runs over all chiral 
multiplets. 

The supercurrent is conserved dfj,S^ = as a consequence of equations of motion. After fixing the 
vector gauge symmetries in the usual way, the vector equations of motion change due to the gauge- 
fixing terms and to the ghost current. The ghosts are scalars under supersymmetry (in particular, 
they do not have superpartners and they couple only to the gauge field but not to the gaugino), and 
one could be worried that SUSY is broken by the gauge choice. Indeed the supercurrent divergence is 
no longer zero, however it is BRST exact (see e.g. [33]). Thus the amplitude for longitudinal gravitino 
emission still vanishes, and the gravitino gauge invariance is preserved. 

The terms proportional to are sometimes omitted from the supercurrent expression (A.28), 
because they do not contribute to the massive gravitino production due to the on-shell condition 
'j^^^ = 0. However, one should be careful to keep these terms if one wants to use the equivalence 
theorem and the massless gravitino gauge invariance, because the supercurrent is no longer conserved 
if they are omitted. 



Goldstino couplings 

According to the equivalence theorem discussed above, the spin ±1/2 component of the massive 
gravitino at high energies can be replaced by the Goldstino coupled to the divergence and trace $ 
of the visible sector supercurrent with coefficients given in (A. 18). The divergence d^S^^^^ measures 
the SUSY breaking in the visible sector, which at energies lower than the messenger scale looks like 
explicit breaking by soft terms. In absence of soft terms d^S'^-^^ = as a consequence of equations 
of motion. Nonzero soft terms modify the equations of motion, so that d^S^-^^ ~ mgoft / 0. For 
dimensional reasons we can neglect dimension 2 soft terms (i.e. scalar squared masses): only soft terms 
with dimension 1 (i.e. gaugino masses M and trilinear scalar couplings A) contribute to Goldstino 
production at dominant order, 7 oc T^. By taking into account how the relevant soft terms modify 
the equations of motion of particles and sparticles we get 

^M^vis = -'-fKp[Y,Yh'>^'' - M<7(rr»A" + V2[{Awum + iAw)*{cp*)^},] (A.29) 

where i runs over all chiral multiplets and a sum is understood over the components of the gauge 
group. The presence of the second gauge term was first noticed in [34], and is here reobtained via a 
simple direct computation. Notice that to get it, it is crucial to keep the last term in eq. (A.28), that 
does not contribute to massive gravitino production due to the on-shell condition 7^^*^ = 0. 

In the MSSM the relevant soft terms are the three gaugino masses Mi^2,3 and the top A-term, At: 

■^soit = Yl ^^N^N + XtAtiQUH^ + h.c.) + • • • . (A.30) 

N=l 

^^See [32], page 141. Notice that [32] has a misprint in normalizing the second line of the RHS of (27.4.40), 
cf (26.7.10). The difference in 7^ in the terms involving gluinos is because our gluino-squark-quark couphng is 
real: Aour = ij^^his- The extra i is then compensated by the difference in 7^^, = i7(^jg. 
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Finally, we elaborate on the Goldstino coupling to $, finding that it can be neglected in the MSSM. 
Using -i^'lYn^h^ = 0, (A.27) implies 

^M^M = _^ [-2(Z?>,)V.ei + ^iW*{(t>*%] - ^ig{<P*^<PjyX'' (A.31) 

Rewriting the first term as 

and using the fermion equation of motion, several terms cancel and we remain with 

a''s^' = 2V2d^i(j)*a^^i) + 2V2i[^*W*j - 2W*i<l)*Mj (A.32) 

The first term does not contribute to massless Goldstino production rate since xdui4'i'^u^i) vanishes 
on-shell due to Oyd^x = 0- The second term vanishes if Wj is a quadratic function of the fields, i.e. 
for cubic terms in W . We thus conclude that the only nontrivial coupling to Goldstino arising from 
the (A. 15) vertex is due to the ^-term and is of the form 

^ x{HiH2 + HiH2) + h.c. 



Mpi 



This vertex is irrelevant at energies much bigger than /i. 

The reason for the above is that the trace of the supercurent $ falls into a supersymmetric "anomaly 
multiplet" 

where is the trace of the energy-momentum tensor expressing the scale invariance of the theory, 
and is the current of the i?-symmetry under which all chiral multiplets have charge 2/3 (see [32]). 
In MSSM, both the scale invariance and the above i?-symmetry are broken classically only by the 
//-term, and this explains $ ^ fi. At quantum level the scale invariance and the i?-symmetry are 
anomalous, e.g. df^R^ is given by the triangle anomaly equation: 

a ^N9% p(N)p(N) 

where the anomaly coefficients = {11, 1, —3} are the same as the one- loop /3-function coefficients 
of the MSSM gauge groups, which is again related to the fact that d^R^ and are in the same 
supermultiplet. Since supersymmetry relates $ to dfj_R^, one can show that (see [35]) 



bNg_ 
levr 

Below we argue that this equation can be used also at finite temperature. 



Y.-^PPb'.l'W^'''^. (A.33) 



A. 4 Gravitino and goldstino couplings at finite temperature 

Gravitino production from a supersymmetric thermal plasma is best studied in terms of its non- 
time ordered propagator 11^ (P) given by eq. (2.2). Supersymmetry is broken by finite temperature, 
but this breaking is spontaneous, so the supercurrent remains conserved: 9^5^ = holds as an 
operator equation. This means that the production rate of longitudinal gravitinos vanishes also at 
finite temperature. Equivalently, 11^ (P) is invariant under gauge transformations of the gravitino 
polarization tensor, dl^^lu = P^Ay + PuB^, with arbitrary A,B. The statements of the previous 
paragraph should hold identically in any computation including all diagrams to a given order in the 
thermal bath coupling g. In practice, however, it may be difficult to see the vanishing of 511^ explicitly. 
E.g. as explained in section 2 we are resumming a well-defined class of physical effects to order g^: 
those enhanced by a 1 ^ 2 phase space factor, unlike a generic 0{g^) correction. More precisely, 
diagram D is computed including thermal corrections to the propagators of particles to which the 



29 



gravitino couples, while diagrams 81^2,3 are computed using tree-level propagators. In particular, we 
do not include corrections to the gravitino vertex. For this reason we expect a residual gravitino gauge 
dependence, which we believe to be of relative order g^/vr^ with respect to our result. The reason is 
that in our calculation of the massless gravitino production rate, thermal masses act similar to soft 
SUSY-breaking terms, modifying equations of motion by terms of order g^T, so that d-S ^ g^T rather 
than being zero. This means that JII"^ ~ 0{g'^). We see that this non-gauge invariance is of the right 
order of magnitude to be cancelled by vertex corrections. The above residual non-gauge invariance 
can be tolerated when computing the massless gravitino production rate. 

When computing the Goldstino production rate, we have taken into account that, in absence of 
soft-SUSY breaking, the Goldstino coupling to d ■ S vanishes at finite temperature by evaluating the 
divergence of the supercurrent before computing the thermal matrix element, i.e. we start the finite- 
temperature computation from eq. (A. 29). Since we do not evaluate vertex corrections, this procedure 
is expected to give a result with the same 0{g'^ /ir'^) error as the gravitino production rate. 

Finally, the anomaly relation (A. 33) valid at zero temperature also holds at finite temperature. The 
argument is the same as in case of the supercurrent conservation: the thermal bath is a background, 
and (A. 33) is a dynamical property of the Hamiltonian valid for any background. In practice this 
means that the Goldstino coupling to $ can be neglected. 



B Vector propagator at finite temperature 

We list the full one-loop expressions for thermal corrections to a vector with four-momentum K = 
{uj,k) (K'^ = uj'^ — k'^) with respect to the rest frame of the thermal plasma. In general, we denote 
by Ufj, the four- velocity f/^ of the plasma. We use the Feynman gauge where all effects are condensed 
in two form factors even in the non-abelian case [9]. Polarizations are conveniently decomposed in 
Transverse (i.e. orthogonal to K and to k), Longitudinal (i.e. orthogonal to K and parallel to k) and 
parallel to K. The corresponding projectors (11-^ -|- 11^ -|- n^)^i, = —1]^,^ are 

uT ~ , ^^k, /o \ . 



K = -^.. + ^-n^., (B.ib) 



where fj^y = r]^i, — U^Ui,-, = — {K ■ U)U^. The vector propagator is 

*D 



— TTq — TTj- — T^O — T^L K"^ 



(B.2) 



In the following ~ denotes the HTL limit, where the result can be expressed in terms of the vector 
thermal mass rriy = ^g'^T'^{N + Ns + Np/2), where the Sector, Fermion and Scalar coefficients 
are defined having in mind a group SU(A^) with Np massless Dirac fermions and Ns scalars plus 
anti-scalars in the fundamental representation. Table 2 lists the explicit values of N, Np, Ns in the 
SM and in the MSSM. The one-loop quantum correction at T = in the MS scheme is 

^^2Nf + Ns-5N -K^ 
^0 = 9K ^^-2 In -3^ (B.3a) 

where the gauge-dependent vector loop gives a negative contribution to spectral densities above the 
light cone, at > 0. The thermal corrections are 
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Standard Model 




MSSM 


Vector 


N Np Ns 


N 


A^^ Ns 


Gluon SU(3)c 


3 6 


3 


9 6 


Weak SU(2)i 


2 6 1/2 


2 


9 7 


Hypercharge U(l)y 


10 1/2 





11 11 



Table 2: Numerical coefficients for vector thermal mass rriy = ^g'^T'^{N + A^^ + Np/2). 



where 



Gs 
Gf 
Gv 
Hs 
Hp 



-^g\NsHs + NfHf + NHy) ^ -^{L + iWv, 
- Y + ^i^sGs + NfGf + NGv) ^ml{l + 



H 



V 



dp 
2^ 
dp 
2^ 
dp 
2^ 
dp 
2^ 
dp 
2^ 
dp 
2^ 



4p L_ 

\ K\ ' 
^ 2k 



rp2 

nsip) ^ — 
o 

nB{p) - 



M k 
2pL + — + -L_ 
k A 



M' 



npip) ^ 



■ M k^ 



12 

nsip) 



having defined lo± = (w it A;)/2, 

L = 1 



-In — , 

k UJ- 



P + UJ- 



M = {p + uj+){p + LO-) h-1 



P + UJj, 



{p — uj+){p — UJ-) In 



P — LO- 

P-UJ+ 



2kp 



B.3b) 
B.3c) 

B.3d) 
B.3e) 
B.3f) 
B.3g) 
B.3h) 
B.3i) 



p + LO- ' ' ' " ' p — OJ- 

See [9, 10] for previous results. We added scalar loops and wrote logarithms such that imaginary parts 
(needed to get spectral densities) are obtained using the prescription a; ^ a; + iO^, with Inz having 
a cut along the negative real axis. We emphasize that our expressions cannot be simplified using 
ln(a6) = In a + In 6, because this would give wrong imaginary parts. The spectral densities employed 
in eq. (3.5b) are defined as 

1 1 
PT = -2Im^r,^— — , p^ = _2Im— — . (B.4) 



K"^ — VTo — TTt ' k'^ K"^ — T^O — T^L 

As well known [9, 5] pL contains a collective longitudinal excitation, that corresponds to longitudinal 
waves of electric fields allowed by Maxwell equations with vanishing dielectric constant. 



C Fermion propagator at finite temperature 

Fermions can receive thermal corrections from gauge and Yukawa couplings. In HTL approximation 
the full result is determined, in a generic non-supersymmetric theory, by one parameter, the thermal 
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mass: 



'Cr 



9' + 



16 



(C.l) 



where we used for g, A, Cr the same notation as in (3.2), except that g, A here denote non-supersymmetric 
couphngs. The same parameter mp controls the full one- loop expression in the Feynman gauge. 
See [8, 10] for previous results. The spectral densities for particles (p+) and holes are given by 



P± 



— Im 



1 



A^ 



In- 



_}_(^ 2 I 

2^2^ 8 ^ 16'"^ /i2 



where the T = contribution gives a spectral density only above the light-cone, and 

Loj^ + Ll!± 



ku> 



(C.2) 



(C.3) 



The functions oj± = {oj ziz k)/2 and L± are the same previously defined for vectors. The last terms 
is the HTL contribution (complex only below the light-cone). Again branch cuts are defined by the 
prescription u> ^ uj + iO+. 
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